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Abstract 



A quantum statistical model of nuclear multifragmentation is proposed. The recur- 
rence equation method used within the canonical ensemble makes the model solvable 
and transparent to physical assumptions and allows to get results without involving 
the Monte Carlo technique. The model exhibits the first order phase transition. 
Quantum statistics effects are clearly seen on the microscopic level of occupation 
numbers but are almost washed out for global thermodynamic variables and the 
averaged observables studied. In the latter case, the recurrence relations for mul- 
tiplicity distributions of both intermediate-mass and all fragments are derived and 
the specific changes in the shape of multiplicity distributions in the narrow region of 
the transition temperature is stressed. The temperature domain favorable to search 
for the HBT effect is noted. 
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I. INTRODUCTION 



A considerable interest is paid to heavy ion collisions which result in the formation of 
residual nuclei with the excitation energy about 5 — 10 AMeV. In this range of excitation 
energies one can probe an intriguing possibility of the nuclear liquid-gas phase transition. 
Below these energies the produced system is deexcited by the evaporation of particles, mainly 
nucleons, as would be expected from a hot liquid drop, whereas at higher excitation ener- 
gies the excited source explodes with simultaneous formation of many Intermediate Mass 
Fragments (IMF). The way from the production of a residual nucleus to the formation of a 
multi-fragment system (the multifragmentation process) is a complex phenomenon where a 
relative role of dynamical and equilibrium effects is still not completely understood. Such a 
situation gives rise to the variety of approaches for describing the many-facet nuclear mul- 
tifragmentation : from an evaporative picture [|1] and simple percolation models till a 
refined statistical treatment f§-[7| and complicated dynamical models [§-412[] (see the cited 
above review-articles and references quoted therein). 

In this work we develop the Quantum Statistical Model (QSM) of multifragmentation 
based on the canonical formalism. Quantum statistical treatment has been applied earlier 
to the multifragmentation process in Refs. fT3|-|T7|| but within the grand canonical ensemble 
providing the conservation of proton and neutron numbers in the mean. The canonical 
method, which is more suitable for the description of multifragmentation of a finite system, 
corresponds to the exact conservation of baryon number. However till now it was used 
only for the case of classical Boltzmann fragments [|6| |T8| , |T9"f . The mathematical basis of our 
consideration is the recursive technique |20| which allows one to eliminate the Monte Carlo 
calculations |J and to present results in an analytical form being rather similar to that 
for the Boltzmann case [ 18|, 19 1 . The model developed in this paper addresses the issue of 



quantum effects in the multifragmentation process. Although these effects a priori cannot be 
large, in some particular cases they may be important. For example, it was demonstrated 
that the yield of Li-isotopes, which is used in the studies of the caloric curve behavior 
to determine the system temperature by means of the isotope-ratio method, is sensitive 



to quantum statistics |]T7] . Intensity interferometry for identical fragments may give 



certain information on the size of fragmenting sources what is of primary importance for 
understanding the multifragmentation process. It is also well-known that Bose-Einstein 
correlations can be the source of salient short-ranged phase-space correlations p2]j23| . In 
general, the study of quantum systems of identical particles with incomplete information 
may give intriguing and counterintuitive results even in simple cases as exemplified in . 

This paper is divided as follows. Thermodynamics of ideal quantum gas of identical par- 
ticles is considered in Sec. 2. In this section we introduce the method of recursive relations. 
The canonical ensemble approach is generalized in Sec. 3 to the case of mixture of differ- 
ent quantum gases for a given composition of its components. Based on these results, the 
QSM is introduced in Sec. 4. Technical details of calculations are given in the Appendixes. 
The partition function, various thermodynamic quantities as well as multiplicity fragment 
distributions and their moments are derived analytically via the recurrence relations. Il- 
lustrative numerical results of the QSM are also presented and discussed in this section. 
Finally, concluding remarks are given in Sec. 5. 
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II. PERFECT GAS OF IDENTICAL PARTICLES 



Let us consider an ideal gas of N identical particles, either bosons or fermions. We 
want to find the canonical ensemble partition function and the ensemble averaging of par- 
ticle multiplicity. Our consideration will be based on constructing the recurrence relations 
which enable an exact calculation of both the statistical partition function and the ensemble 
averaging [2Cfl. 

Let a system of ideal bosons (fermions) be enclosed in the volume V at the temperature 
T. Without any loss of generality we can limit ourselves to the case of spinless particles, 
i.e., bosons and fermions differ only by the symmetry of the wave function with respect to 
permutation of identical particles. A microscopic state of the system is determined by a 
set of occupation numbers of single-particle levels {vp\, so in the considered state there are 
vp particles with momentum p | 2~5| . Single-particle energies are given as ep = p 1 /2m where 
p(p( x \ p( y \ p( z ') is the eigenvalue of the momentum operator of a single particle 



,(«) 



pv», = A p W n (a) ; n w = o, ±1, ±2, . . . , a = x, y. : 

and Ap( a ) is a step of the momentum grid towards the direction a 

2nh 



,(«) 



Ap 



(«) 



(1) 



(2) 



In the second quantization method, Eq. ([J) determines an infinite set of discrete values 
of momentum components of p. Every discrete value of p corresponds to a single-particle 
state with energy ep and to a particle number up. The total particle number N and the total 
energy of this microscopic state E^^yjy are 



N = J2 vp 

v 

E {isp}N = T. e v v v 

V 

The allowed occupation numbers are 



^=0,1,2, 
up =0,1 



N (for bosons) 

(for fermions) . 



(3) 
(4) 

(5) 



The canonical ensemble partition function for the perfect Bose (Fermi) gas of iV identical 



particles is 25 



Qn = E e 



(6) 



where the total energy Ef n | ^ is defined by Eq. the Lagrange multiplier /3(= l/T) and 

the sum runs over all sets of occupation numbers (0), taking into account the conservation 
of the total particle number N. The expectation value of an operator O over an ensemble 
of iV identical particles can be found in a standard way: 



3 



Q 



<o>„ = ^- E o w «^W". (7) 



The partition function (^) can be treated as a generating function for moments of occupa- 
tion number distributions. Then the average value (|7|) for different moments of occupation 
number distributions can be found in a following way. Let us introduce a new variable 

xp = e'fcp (8) 

in the partition function @ 

qn= e n • ( 9 ) 

V 

Then, the mean value and highest moments of occupation number distribution can be ex- 
pressed through the derivatives of the partition function (||) with respect to xp 



1 d d 

Q N v dxp v dxp> 
1 d d d 



{vp Uf) N = — xp — xp ^— Q N (11) 



^ V * ^' )n = Q~ N ^ dx~p X * -dx~p x ? Qn ■ (12) 

The frequently used dispersion of the occupation numbers is 

A(up up) N = (up up) N — (up) N (i>p)n 
d d 

= xp — xp - — In Qn . (13) 

OXp OXp 

Let us show now that the partition function Qn in the occupation number representation 
(§) is reduced to the multinomial sum [^,^] which can be solved by means of recursive 
equations. The comparative analysis of explicit expressions for the partition function Qn 
and the mean multiplicity (vp)N for two systems having N and (N — 1) particles reduces to 
the following relations (see Appendix A) 

( u f)j Qj = e~@ £ P Qj-i [I =F (vf)j-i] , ] = l,2,...,N (14) 

with the initial values : Q = 1 and (vp)o = 0. The sign " + "/" — " in Eq. flip) corresponds 
to the Bose-Einstein/Fermi-Dirac statistics, respectively. Eq. flT4] ) is valid for both finite and 
infinite set of single-particle levels. Let us multiply Eq. (14) by : (=Fl) Ar_J ' e -P £ p(N — j) ^ 



and sum both left hand side (l.h.s) and right hand side (r.h.s) of this equation over j from 
1 to N. This results in the following representation for the mean occupation numbers : 
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iV 



j=i 



(15) 



By the substitution of (|lq) into the total particle number conservation condition (y), we 
arrive at the related set of recurrence equations for the partition functions (0) of systems 
with different number of particles 



N 



Q 



N 



N 



E fj Qn-j 



Qo = i 



J'=l 



where 



(16) 



(17) 



This solves the problem of finding the canonical ensemble partition function. The sum in ( |T7D 
runs over the whole infinite set of eigenvalues of the momentum operator of a single particle. 



Inserting different Q?' s m to the recurrence relations fllq) gives rise to the multinomial sum 



Q 



N 



E 

N 

{rij}, E 

.7 = 1 



N 

n 



(18) 



iV 



which is the well-known quantum-group Meyer expansion for the ideal Bose (Fermi) gas of 



identical particles in the canonical ensemble method (rij are integer numbers) 



One 



should note that the recurrence relations (|I6"D are satisfied for any multinomial sum (WE 



independently of an explicit form of the fi variable p0| . The recurrence relations for the 
case of the ideal gas of N identical fermions were obtained recently by Pratt [57] using of 
somewhat different technique. 

The mean occupation number is given by Eq. (|T5|). To get the second moment of the 
occupation number distributions, let us compare the r.h.s. of Eqs. fll5|) and ([[(]). Using the 
definition @ of x$ variable, we find the following relation between the first derivative of the 
partition function Qn with respect to xp and the partition functions of systems with smaller 
number of particles 



dQ 



N 



N 



8Xr 



E (Tiy +1 4 Q 



N-j 



(19) 



Substituting 



°p i=i 

|^) into Eq. (|TT|) and using again (|J), we obtain finally 

-. TV N—i 

<W*)n = 7T E E (Tiy + > e-P** e-M 1 * Q N -i-j 

V7V i-x j=1 



+ <W TT E {Tl) j+X e-^jQ 



N-j 



(20) 
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The knowledge of the mean occupation numbers (|15D, their second moments (|20| ) or disper- 
sions flT3] ) is sufficient for evaluating most of the thermodynamic quantities. Higher moments 
of occupation numbers can be found, if necessary, by the same procedure as presented above. 

Similar results can be obtained for the limiting case of an ideal gas of classical Boltzmann 
particles in the canonical ensemble. One should take into account the "correct Boltzmann 
counting" |p5| , i.e., that for an arbitrary set of occupations numbers satisfying (|3D there 
are generally more ways to construct an ensemble as compared to the Bose gas because the 
interchange of any distinguishable particles leaves the set {up} unchanged 



Qn 



1 

iV! 



E 



N 



M 



-/3E 



{u f }N 



(21) 



and the sum runs over all allowed values 

Lto account the conservation of the total 
This sum is the multinomial Newton sum which can be calculated 



where the total energy Ef v 1 jy is given by Eq. 

of occupation numbers : up — 0, 1, . . . , N, taking into account the conservation of the total 
particle number 
exactly [[25 



Q 



l 

M 



N 



E 



-pep 



(22) 



The partition function for the ideal Boltzmann gas (|22|) corresponds to a single term of the 
multinomial sum (IS) of quantum Bose gas with n\ = N and rij = for j ^ 1 



N 



Qn ~1v\ 



E 



-p£p 



(23) 



Therefore, all results presented above for the Bose gas and expressed through functions fj 

and e~ftj £ P, can be immediately applied to the ideal Boltzmann gas after a simple substi- 
tution 




5j i uj 



J = 1,2, 



N 



(24) 



where 8j t i is the Kronecker symbol. Thus, in the limiting case of the classical Boltzmann 
statistics only the term with j = 1 survives in each sum over j in Eqs. ([15])-(|17]), (|19D , (^0|). 
For example, inserting ( |24]) into Eqs. (|16|) results in the following recurrence equations for 
an ideal Boltzmann gas 



Qn — j~f u Qn-i 



from where directly follows Eq. 
occupation number 



Using flU), (H) and 



N-l 



Qn 



N 



-Ph 



E e 
v 



-pEf 



(25) 

we have for the mean 
(26) 
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Analogously, the second momentum of the occupation number distribution reduces to 
the form 



= N(N - 1) 



-(3e$ e -f3£f 



J2 e~@ £ P 



~2 +5pf N 



E e 

v 



(27) 



which is well-known in statistical mechanics p5[ 



Let us now consider fluctuation of particle number distributions over different quantum 
states. The mean squared occupation numbers (f|)jv for Bose and Fermi statistics are given 



in Eq. fl20j) which can be rewritten as follows 



1 JV a- 

rT E (Tl) J+1 e~^ e ? Q N ., [T(j - 1) + j] ■ 



(28) 



Using Eqs. (fig ) and (p8|) we get for the occupation number dispersion in the Fermi gas 

((Az/p-) 2 )^ = ({vp- (v?) N ) 2 ) N = (v f ) N [1 - (vjr) N ] . (29) 

This expression coincides exactly with the appropriate expression calculated within the 
grand canonical method p8| . This is due to the equality : (^|)jv = (^p)n, which follows 
from (|j~5l) and (p8p. 

Similar consideration for the Bose gas gives 



JV 

O £ 



<(A^)% = J- £ ( T l)^ +1 e -# e ? (2j - 1) Qiv-i - (u p 



(30) 



For the case of the Boltzmann gas we have from Eqs. ( ffif) and (27) 



<(A^) 2 >£ = (u f ) 



B 
Pi N 



Mi 

iV 



(31) 



Here the first term corresponds to the grand canonical result |28[ : ((Au^) 2 ) GC = (Pp)GCi 
where (■ ■ -)gc denotes the grand canonical ensemble averaging. As is seen, in the canonical 
ensemble treatment an additional negative term arises which vanishes in the thermodynamic 
limit iV — > oo. 

Let us introduce the quantity 77 as an integral measure of correlations between the occu- 
pation numbers p8 



p 



E (^) 



(32) 



In application to the Fermi gas, using Eq. (p9|), we have 



-]y E 



(33) 
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It is seen that rj is negative for all values of T, V and 77 = —1 at the limiting temperature T = 

0. In the grand canonical ensemble for the Fermi gas we have : rj = — (1/ (N)ac) XX^gc ■ 

v 

For the correlation of identical particles in the Bose gas, substituting fl3~0|) into (|32|) and 
using (P~5|) , we have 

V = ^ 1^ E fi 2(i - 1) Qtf-i - E • (34) 

At T = one finds : rj = — 1, because : Qtv = 1 and /j = 1 for any N and j . In the grand 

canonical ensemble : 77 = +(1/ (N)ac) { v p)'gc • Similarly, in the Boltzmann limit of the 

p 

Bose gas with using Eq. (|32|), we have : 

V B = ~jp E ("*)% ■ (35) 
p 

In this case the grand canonical method yields rj = 0. At T = one finds : rj = — 1, because 
the mean occupation numbers are equal to : (z/p*) n = NSp^ . 

Thus, the occupation number dispersion in the canonical ensemble for Bose and Boltz- 
mann statistics differs from those in the grand canonical, because the statistical indepen- 
dence of a particle in a state p is broken due to the requirement of the total particle number 
conservation (Q). At the same time, all three statistics predict for the value of the correlation 
coefficient 77 — > — 1 when T — > 0. 



III. IDEAL GAS MIXTURE OF BOSE-FERMI FRAGMENTS 



Let us consider a system consisting of Aq individual nucleons, N 2 clusters (fragments) 
made of two nucleons and so on until Na fragments of A nucleons, contained in the volume 
V at temperature T. The composition of this mixture is characterized by a given set of 
fragment numbers (Aq, N 2 , . . . , Na) and satisfies the nucleon number conservation condition 

A 

i N = A . (36) 

i=i 

A fragment of any species i follows either Bose-Einstein or Fermi-Dirac statistics and a 
number of fragments of a given species N is assumed to be conserved at all permutations. 
Let us derive the thermodynamic relations for such a mixture of ideal Bose and Fermi 
fragments within the canonical ensemble method. 

A microscopic state of an ideal Bose or Fermi gas subsystem of i-th species is defined by 
a set of occupation numbers {vpi\ for single-particle levels of the fragment i. This implies 
that the microscopic state in question has clusters of mass % with momentum p ||25| . The 
single-particle energies are : epi = pi 2 /2rrii. 

The fragment number A^ of every species % and the total energy of quantum gas mixture 
of a given partition or composition : E^^yjy^ ^ jy^, are defined by the following 

expressions (cp. to @ and @ ) : 
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Ni = J2 upi , i = l, 2,..., A (37) 
p 

E {v&}N 1 ,N 2 ,...,N A = T i £ pi v Pi (38) 

i,p 



with the allowed values of occupation numbers given in fl5|). 

Generalization of the canonical ensemble partition function (|6|) to the case of an ideal gas 
mixture of non-identical Bose- Fermi particles with the conservation of a number of fragments 
of every species ([37]) can be written as [p5] 

W A (N 1} ...,N A )= E e-^^u-.NA (39) 
where the sum runs over all partitions of occupation numbers (|5]) 

E = E n 5 (e (40) 

satisfying the conservation laws (p7|) . Generalization of ensemble averaging (^) to the case 
of an ideal gas of Bose-Fermi fragments of different masses is 

(0)n u -,n a = w (N l N , E 0{N 1} ...,N A ) e - pE MNi,... t N A (41) 

Since fragments do not interact with each other, the total partition function (|39| ) is reduced 
to the product of partition functions Q(Ni) for subsystems of identical fragments 

A 

W A (N 1 ,...,N A ) = l[ Q{Ni) . (42) 

i=l 

Evaluation of the partition function Q(Ni) for a subsystem of iVj identical Bose (Fermi) 
fragments of the mass i has been considered in the previous section. Then, Eqs. @, ([16]) 
and ( |T7| ) can be rewritten in the new notation as follows 

Q(Ni)= E l\_e-P £ & u & (43) 

V 

= ? E fuQ(Ni-l) , i = l,2,...,4 (44) 

/« = E (45) 

p 

where the signs (+) and (— ) correspond to the Bose-Einstein and Fermi-Dirac statistics, 
respectively. 

In the same manner, the mean occupation numbers and their second momenta for a 
subsystem of iVj ideal Bose (Fermi) identical fragments of the species % can be expressed by 
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the canonical ensemble partition function ([43|) for smaller number of fragments of the same 
species (see ©, (H) ) 



= E Q(Ni - I) (46) 



fa* = 7vm S E (Tl) /+fc e"^ e"^ Q(iV, -l — k) 

Vl iv «J ;=i fe=i 



+ V g^y E (Tl)' +1 ZQ(iV 4 - . (47) 



By the subsequent substitution of the recurrence relations ( f4"4|) into Eq. (f42|), the parti- 
tion of an ideal quantum gas mixture is reduced to the generalized quantum-group Meyer 



expansion (cp. to Eq. (TLB 



W A (N U ...,N A ) = n E n ^TT • (48) 

Comparing the definitions of ensemble averaging for identical particles (0) to that for 
a mixture of Bose- Fermi fragments (f4l|) , we can relate the first and second moments of 
occupation numbers as well as their dispersions, to appropriate quantities for subsystems of 
identical fragments: 

(v&)N lt ...,N A = {Vpi)N x (49) 
{vpi Vfj)Ni,...,N A = Sij {Vpi ^pij)Ni + (1 — Sij) {vtfijNi {vpj)Nj (50) 
A(z/p-j Vpj)Ni,...,N A = ( v ft u p'j)N 1 ,...,N A — (vpi)Ni,...,N A ( u p'j)N 1 ,...,N A • (51) 

Therefore, to find the moments of occupation number distributions (f4"9"|)-(|5"I]) for a mixture 
of quantum gas with fixed composition of fragments, it is sufficient to know the occupa- 
tion number moments for subsystems of iVj identical fragments. The latter are related via 
Eqs. (|46|), (f47|) to the subsystem partition function Q(Ni) and can be calculated recursively 
by means of Eq. ([44]) . 

Let us turn now to the calculation of general quantities characterizing thermodynamic 
behavior of a system. By definition, the average energy of an ideal gas mixture of Bose- Fermi 
non-identical fragments is a partial derivative of the logarithm of the partition function 
with respect to the temperature at a constant volume V 



E(N l ,...,N A )=T* d -± W ^-' N ^ 



dT 



( E {N j z},N h ...,N A )»i>-,NA (52) 



where the total energy expectation value is calculated over the ensemble ([y]) and the sum 
runs over all values of the occupation numbers (§). Because the gas is ideal, the average 
energy of the system is the sum of subsystem energies 

E(N h ...,N A ) = J2 £ & Mn u ...,n a ■ (53) 
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Inserting the mean occupation numbers fl49|), expressed through the partition functions for 
subsystems of identical Bose- Fermi fragments (|46|) , one obtains 



E(N 1 ,... 1 N A )=J2 



Ni 



E Wi Q(Ni - I) 



(54) 
(55) 



The summation procedure over the fragment momentum p in Eq. ( p5|) is presented in the 
Appendix B. 

By the definition, the pressure in the system is 

BE, 



P{N 1 ,...,N A ) 



d\nW A (N 1} ...,N A ) 



dV 



J {N fl },N 1 ,N 2 ,...,N A 



Inserting here the average energy ([38|), we find 

P(N 1 ,...,N A ) = 



Be- 

E p* / \ 
\ V Pi)Ni,.. 



%,p 



BV 



BV 



,N A 



,N A 



(56) 



(57) 



Making use of Eqs. (fig), (|49|), the pressure in a mixed system (|57j ) can be related to the 
partition functions for the Bose-Fermi subsystems of identical fragments ([44]) 



P(iV 1; ...,iV A )=E 



E p ^ W-O 



^ = -fri) m E ^ e_/3 ^ 



(58) 
(59) 



In a similar way, one can calculate the specific heat of the system which is defined as a 
derivative of the average energy (p2]) with respect to temperature at the constant volume 



C V (N U ...,N A ) 



8E{N U ...,N A ) 



dT 



v 



N 



a )n u ...,n a - ( e {v pi },N 1 ,...,N a )nu..,n a 



J^2 



(60) 



Inserting (|38|) into ([HP and using (^9j)-(|51f) , we express the specific heat for mixture of 
an ideal quantum gas as 



Cy(Ni, . . . , N A ) = p 2 E ^{NpiVpj)Ni,...,N A £pi £fj 



(61) 



Hence, the specific heat can be expressed via the partition functions of identical Bose-Fermi 
fragment subsystems (|4"5D, if one applies Eqs. fl4"6]), ( f4"T| ) and (|^]), (pH|) 



c v (n 1 ,...,n a ) = p 2 '£ 



E £ * 



I £ik 



Q{Ni -l-k) 



Q(Ni) 



Q{Ni - l)Q(Ni - k) 



+ P E 



i 



(62) 
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In the above equations, the quantity en is defined by Eq. (|55D and the average squared 
energy of a fragment is 



?. = rzpi)^ 1 £ 4 e -/5fe pl (63) 



Evaluation of average values of e~a and e\ is given in the Appendix B. 

Let us consider now the Boltzmann limit of an ideal gas mixture of Bose- Fermi fragments. 
The partition function of the ideal Boltzmann gas of non-identical fragments taking into 
account the conservation law for fragment numbers of every species (^) is (cp. to (^)) : 

W° A (N U ...,N A ) = ^— £ nf W^-^},^, Wj (64) 

n N k \ {„„} *-i VV "*7 

fc=l 

where total energy -^{^.j _/Vi _/y 2 _/y A is given again by Eq. (j38| ) and occupation num- 
bers take values : = 0, 1, . . . , iVj. Similarly to the quantum partition function (|39D , the 
classical partition function (B3) can be expressed as a product of partition functions Q B (Ni) 
for subsystems of iVj identical Boltzmann fragments of mass i. As was shown in the pre- 
ceding Section, the partition function Q B (Ni) for the subsystem of iVj Boltzmann identical 
fragments of mass i, which is obtained from the appropriate quantum equation by the simple 
substitution of (|24|) 

Q B (N t ) = , = fa = E e"^ (65) 



corresponds to a single term of the multinomial sum (|18D with rin = iVj and = for 
j 1 (i = 1, 2, . . . , A) |p5| . So, the substitution of (|24l) into the recurrence relations (fUD 
leaves only one term 

Q B (iV 4 ) = ^^Q B (iV J -l) , z=l,2,...,A (66) 

and the partition function fliSD for the mixed Boltzmann gas of fragments takes very simple 
form: 

A M 

W*(N 1 ,...,N A ) = l[ ^_ . (67 ) 

i=l l ' 

Eqs. (|49|)-(|5TD for the moments of the occupation numbers remains unchanged but the 
evaluation of the corresponding momentum components is simplified noticeably (see (pB| ), 
(1^)). Taking into account Eqs. (|55|) and (|6~5D for the average energy (|54D , we have 

E B (iV 1 ,...,iV A )=E Si^ (68) 

i=l 

where 
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E eae"^ 



p 

From Eqs. (|58|), fl59f ) and fl65f ), one finds that the pressure in a mixture of the ideal Boltzmann 
gas is 

A 

P B {N 1 ,...,N A ) =E P'iNi (70) 

4 = 1 

where 

P. = - '' 3V ■ < 71 > 

p 

Specific heat of the system in the Boltzmann limit is obtained from Eqs. (j62|) , (|63|), ( E5[ ) and 



^(iV 1; . . . , iV A ) = /5 2 E ( £ ? - 5" 2 ) iV, (72) 

i=i 

with 



12 _ p 



P 



(73) 



These relations will be used below when describing the thermodynamics of the quantum 
statistical model of the multifragmentation. 



IV. QUANTUM STATISTICAL MODEL OF NUCLEAR 
MULTIFRAGMENTATION 

A. Model formulation 

Based on the recurrence technique discussed above, we shall formulate now a quantum 
model of fragmentation of a finite nuclear system into nucleons and nucleon clusters. This 
quantum fragmentation model is constructed in the the framework of the canonical ensemble 
description in the occupation number representation. Along this way, we shall closely follow 
the physical approximations used in the popular statistical fragmentation models in the 



canonical ensemble |q, |29|j30|| . In particular, the partition function of a fragmenting system 
is defined as a sum over all fragment partitioning, each terms of which is taken with the 
weight associated with the partition function of a given fragment composition. 
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The problem of the partitioning of an integer A into a sum of integer numbers is known 



in the number theory for a long time [P0| , Pq| . Some multifragmentation models pl|j32 
even employed this mathematical result for the description of mass fragment distributions 
assuming equal probability for any fragment composition. A partition of A nucleons into 
clusters (N\, . . . ,Na) corresponds to the fragmentation of a nucleus A into Ni individual 
nucleons, N2 fragments made of 2 nucleons etc. with the nucleon conservation constraint 
(|56j). To construct a multifragmentation model, it is necessary to define a probability of 
any partition (iV\, . . . , Na) with the natural condition that the sum of partial probabilities 
over all final states equals one. In the statistical multifragmentation models |J, the non- 
normalized probability of a fragment configuration is given by its partition function. This 
is just the quantity Wa(N±, N 2 , . . . , Na) studied in Sec. 3. 

Strictly speaking, the ideal gas approximation does not address the problem of the nu- 
clear multifragmentation process properly, especially with respect to a possible relation of 
multifragmentation to the liquid-gas phase transition. The physics content of the model can 
be however improved by noting that the decisive condition allowing to develop the general 
scheme of recurrence equations is the additivity of the total system energy with respect to 
subsystem contributions. So, the interaction of i nucleons bound in a single cluster can be 
taken into account by adding the fragment binding energy 6, to the single particle energy: 

^ = |Ue. ■ P4) 

In principle, some mean-field interaction energy depending only on i can be also included into 
Eq. ([74]). In particular, such a form has the Coulomb interaction energy between fragments 
in the Wigner-Seitz approximation 



1 A 2 p 2 
fiW-S _ 71 e 



1 



5 4R \A 



1 



which is used frequently in the statistical description of multifragmentation 0,^,0. 

Repulsive residual nuclear interaction between fragments can be approximately ac- 
counted for in equilibrium statistical mechanics by the excluded volume method p5| , |33| . 
Since a real fragment occupies some finite volume Vi, therefore only a reduced volume is 
reachable for their free motion: 

A 

V f = V-Y,v i N i . (76) 
i=i 

It is then assumed that the system may be still treated as consisting of idealized point-like 
fragments in Vf. It is worth noting that a simple replacement in thermodynamic equations : 
V — > Vf, needs certain care. Appearance of an additional N dependence of the partition 
function may result in the loss of thermodynamic self-consistency. Due to the lack of con- 
sistency some basic relations of thermodynamics, such as those relating a derivative of the 



thermodynamic potential to the particle number or the pressure, may be violated |34]-p6 
Moreover, as suggested by ([TDD, the multiplicity dependence influences the probability of 
a given fragment partition which is proportional to the (micro) canonical weight of a given 
partition and therefore should be properly accounted for. Recent discussions of this problem 
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can be found in Refs. f37H40[. In the QSM, the excluded volume effect is ignored and V is 
treated free variable parameter. 

Thus, the canonical ensemble partition function in the QSM of multifragmentation can 
be written as follows : 

Z A = E W A {N U ...,N A ) 

3 

= E E e -P E {v pi }N 1 ,...,N A (77) 

{N 3 },J2jN 3 =A {vpi} 

3 

where the sum runs over all possible fragment partitioning {N\, N2, • • • , Na) subject to the 
constraint (|36|) . The partition function for a given fragment configuration is described by 
Eq. together with the conservation laws fl57p for each fragment i. Because the partition 
function Wa(Ni, . . . , Na) can be expressed through the multinomial sum (|48|) depending on 
variables (fl5]), we have 



za= h n e n^sr- (™) 

iV,}.V.yV, 1 »=1 {n fe },E fc ^ =Ar » ' =1 

The expectation value of a given operator O over full ensemble of a system decaying into 
ideal Bose and Fermi fragments is defined in the QSM as follows (cp. to flUD) : 

(o)a = 4- E (ok.^ WaC^i,...,^) 



i- £ £ 0{i»} e ^{-»}M,...,^ . ( 79 ) 



Z A 



Let us find the recurrence relations for Za- According to (|79f) , the average occupation 
numbers can be rewritten in the following form 

1 A 
Ma= tt E Mm,..,N A II Q( N i) > i = l,2,...,^ (80) 

where Eq. (fi^) was used for Wa(^Vi, . . . , Na)- The insertion of (|49|) and then ( f46|) into 
Eq. ( p0|) and the exchange of the summation order over I and {Nj} followed by the extension 
of the sum over I to the maximal value N by means of the step- function 8(N — I), gives 

1 [A/i] m A 

Ma = — E E m - I) Q(Ni — OH ■ (81) 

The step-function removes here those terms in iVj for which iVj < /. So, the mean occupation 
number is 
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[A/i] „, [A/i 

pi) A = 



d/*>A = -J-E (Tl) m e-^ E E i|E^-^l Q(Ni-l) YlQW) . (82) 



A 1=1 {Nj^} Ni=l \j=l J jjti 

Changing the summation indices : iVj — > JVj + Z, and noting that Aj = for all j > A — il 
as follows from the Kronecker symbol, we finally have 



i [A/i] A-il 

Ma = jr E e"^ E II QM) 



i [A/i] 

= 4" E Z A _« (83) 

^A , =1 

where Eqs. (|42|), ( |77D were used to get the second equality in (|8~3D . After summing over 
momentum p, the mean number of fragments of mass i in the system of size A can be 
related to the partition functions for systems of smaller size : 

i [A/i] 

(N l ) A = — Yl h 2 A-ii (84) 

A A l=1 



where the variables fu were defined in Eq. fl45|) . Note that 'il' denotes a single index corre- 
sponding to the product of 'i' by '1'. 

Using Eq. (84) the recurrence equations for the partition function Z A can be easily 
derived now from the conservation of the total nucleon number ( |36| ) 

1 A [A/i] 

Za = -t E E i fu Z a-u ( 85 ) 
A i=l 1=1 

where Zq = 1 according to ([f7]). The partition function for an arbitrary A is obtained by 
an iterative solution of the recurrence equations fl35l) starting from Zq = 1 

A [A/i] ( f \n u 

za= e n n • ( 86 ) 

{n u },J2iln u =A <=1 1=1 

The second moment of occupation numbers can be similarly derived from Eq. (|79|) with 
subsequent using of Eqs. (|56|) and (f47|) 



)a = 4~ E frl)** e-^e-^i Z A _ lk _ 3l 



Z A K i 

+ 4~ E e~P l£ & I Z A „ U . (87) 

A i 

Summing over all microscopic states of a fragment gives the second moment for the fragment 
number 

! [A/i] 1 
(Ni Nj) A = $ij tt- E ^ fu Za~u + tt- E f ik fi l z A~ik-ji (88) 

Z A l=1 A fc . 



16 



where all functions Z\, Z 2 , . . . , Z A are calculated from Eq. (|^) with fa defined for the system 
of A nucleons. The knowledge of the first ( [84|) and the second (|88"D multiplicity moments 
allows us to find the multiplicity dispersion 

A(N, Nj) A = (Ni N,) A - (Ni) A (Nj) A . (89) 

Similarly to (29), one can define the dispersion of occupation numbers in a given quantum 
state p 

((A^) 2 ) A = (fa - {^) A f) A = {»l) A - ■ (90) 

The mean squared term (v^) A can be obtained for all statistics from Eq. fl3T|), which reduces 
to 



(4) a = y E (Tiy +1 e-to £ * Z A ^ [ T (j - 1) + j) (91) 

where the mean occupation numbers are defined by Eq. fl8"3"|). So, the dispersion ( p0| ) for the 
case of Fermi gas fragments is 

((Az/ p -0 2 >a = (upt) A [1 - (vpi) A ] (92) 

and coincides exactly with the result for the identical fermions (|29|). For the gas of Bose 
fragments, we have correspondingly (cp. (|34]) ) : 



<(A^) 2 )a = E (2j - 1) Z A _ i3 - (^)i . (93) 

3=1 



Note that the calculation of both the partition function Z A and ensemble averaging is 
carried out with the help of recurrence equations (j55p. However, these equations have still 
other representation related to the recurrence relations (jT6| ) for a system of N identical 
bosons and fermions. By simple algebraic transformations, Eq. (J85Q can be reduced to the 
following form : 

1 — 

Z A = ~ A E °i ^A-i (94) 

A i=l 

where 

°< = E E - &0 fc /w • (95) 

k=l 1=1 

It is seen that Eq. fl9"4|) formally coincides with the recurrence relations (0), for which the 
multinomial sum (18) holds. Therefore, the partition function Z A of the canonical ensemble 
in the QSM of multifragmentation can be presented as the multinomial sum ([Tj^) with the 
redefined variables a« and an evident change : iV — > A. The recurrence equations (|9~4l) can 
be considered as a set of A linear equations with A unknown values Z±, Z%, . . . , Z A : 
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3-1 

jZj — ^2 Oj Zj_i = cij 

8=1 



J = 1,2, 



A 



(96) 



Solving (^) by the Kramer method, we get the solution for the partition function (77) as a 
matrix determinant formed by coefficients 



1 

~A\ 



CLA — — (12 

cla-i A — 1 —cii 
a A _ 2 A - 2 











-a>A-i 

' a A-2 
- a A-3 

1 



(97) 



The partition function Qn for the system of N Bose (Fermi) identical particles (0) can be 
also presented in the form of matrix determinant ( |9"7| ) with coefficients : Oj = fa, and the 
evident change : A — > N, Za — > Qn- Thus, the multinomial sum (|T8|) can be unambiguously 
reduced to the matrix determinant (|9~7l). 

The limiting case of the Boltzmann statistics can also be obtained in the way presented 
above. The explicit form of the canonical ensemble partition function for multifragmentation 
of the Boltzmann fragments is 



L A 



1 

~A\ 



E 



A\ 



E 

- a n N k \ {^} 



k=l 




N, 



n 

V 



f 3E {u r!i }N l ,N 2 ,...,N A 



(98) 



where the total energy E{ Up .} Nl ^ ^ NA is given by Eq. ( pgp and the occupation numbers take 
values upi = 0,1, ... , A^. The ensemble averaging is also defined up to a permutation factor. 



Making use of the substitution (gj) for Eqs. (pq), (p6|), the partition function (|98|) results 
in the multinomial sum for which the following recurrence relations are fulfilled [EUJ 



7B 



A 

e n 



.N, 



NA 



1 

A 



E i U i Z 



B 

A-i 



(99) 



&i = fn in (|99| ) denotes the partition function for ideal fragments of the z-th species (also 
the integer numbers n fl were denoted as iVj). 

The recurrence equations of this type were first applied to the multifragmentation process 
by Chase and Mekjian |18| , |19| . These equations were motivated by the Boltzmann treatment 
of fragments with phenomenologically postulated thermodynamic variable uji (the so-called 
tuning parameter). The partition function Z^ in the form ([}?]) remains unchanged but now : 
aj = iuji, for i = 1, . . . , A. The Boltzmann substitution ( p4|) into expressions for the mean 
multiplicity of occupation numbers (|83|) and fragments numbers (183) gives rise to 



-/3e 



pi 



L A-i 
yB 
L A 



Z B A 



1,2,..., A 



(100) 
(101) 
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For the second moments of the occupation numbers (p7|), their dispersion (0) and the 
fragment numbers (|88D , we have similarly 

(»* = e-P** e-foi %^ + 5 i3 5^ (102) 

^A 



((A^) 2 )a = Ma [1 - Ma] + e~ 2 ^ -4^ (103) 

J A 

(Ni Nj) B A = %P + 6^ 5 p7 u t %* . (104) 

L A L A 

Let us now illustrate the influence of quantum statistics on various characteristics of 
nuclear fragmentation in terms of the QSM. Let us consider three versions of the model 
when all produced fragments are treated as either fermions, bosons or classical Boltzmann 
particles. Note that all fragments are assumed to be spinless, i.e., the spin degeneration 
factor is put equal one. In the results presented here, the Coulomb interaction potential 
is taken into account in the approximation (|75|) and experimental values of the binding 
energies 6j are used everywhere unless opposite is said. Internal fragment excitation is 
neglected. Volume of the fragmenting nuclear system V is regarded as a free parameter 
and its value is given in terms of the volume of a system at the normal nuclear density : 
V/V = Vpo/A. 

The momentum dependence of the mean occupation numbers {i>j%) is given in Fig. 1 for 
different fragments produced in the multifragmentation of system with A = 200 nucleons. 
Bose statistics enhances noticeably the yield of light fragments with low momentum (p < 
150 MeV/c) but this enhancement practically vanishes for IMF's. One should note that 
to get the observable momentum spectra for fragments, the given values of (vpi) should be 
multiplied by the phase factor p 2 . 

Temperature dependence of the mean fragment multiplicity (iVj) for the same system is 
presented in Fig. 2. Though the effect of quantum statistics slightly increases with temper- 
ature T, it is seen in fact only for nucleons. This is caused by the fact that one needs many 
identical particles to get a sizable quantum effect but the mean multiplicity for fragments is 
typically small as can be seen from Fig. 2. Note that the influence of quantum statistics on 
the mean fragment multiplicity is considerably less than that of uncertainty in the choice of 
the volume V of the fragmenting system. The mean multiplicity of the lightest fragments 
grows gradually with increasing violence of the collision, while the IMF's exhibit clearly the 
"rise and fall" behavior observed experimentally [ffl]|. The maximum in the T-dependence 



yB 

A 



of {Ni) of IMF's is slightly shifted towards higher temperatures with increasing size of the 
fragment. 

The "rise and fall" behavior noted above is manifested again in the T-dependence of the 
total mean IMF multiplicity defined as : (M) = EZ G ( N i) (see Fig. 3). All three model 
versions give practically the same results and demonstrate rather strong dependence of the 
maximum on the volume of the system. Due to large contribution of nucleons, the mean 
multiplicity of all fragments : (m) = Z^=i(iVj), is an increasing function of T with abrupt 
growth above the threshold of IMF production. This peculiar structure of (m)(T)-curve is 
manifested stronger for heavier systems, as follows from the same Fig. 3. 

Fluctuations in the region of phase transition behave in a special manner. To approach 
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this issue, one can generalize the correlation characteristic rj fl32|) to the case of fragment 
mixture 



E [ ((A^) 2 ) A - ( Vpi ) A ] 
Vi = I x~~/ \ • ( 105 ) 

E \Vpi)A 
P 

Using fl92|), we obtain for the Fermi gas 

vt = -jj^2 ? ( ^ (106) 

i.e., rji is negative for all values of T, V and 77^ = —1 in the limit T — > 0. Using (p3|) , one 
obtains for the Bose gas 




" (Ni) A 

In the Boltzmann limit, this correlation characteristics is 



^ = Wu I fl2 ^zf ~ ^ {u * ?A } (108) 



f u 2(l-l)Z A - il - y £ Ma\ ■ (107) 



p 



and differs from the Fermi g ( |106|) by a positive correction term. 

In Fig. 4 the temperature dependence of rji is presented for various fragments. Indeed, 
there is a striking difference between the cases of Bose and Fermi statistics and this difference 
is maximal at the temperature corresponding to the maximal yield of this fragment i going 
to zero. As noted above, i] < for the Fermi gas. Beside the lightest fragment i — 1, 
the Boltzmann limit follows very closely the Fermi case and the difference between them 
becomes negligible for heaviest IMF's. Because the nucleon component is dominant and (Ni) 
is rather large, its ^-behavior looks like that for (Nx) identical particles with correlations 
which are positive, negative or zero for bosons, fermions or Boltzmann particles, respectively. 
Deviation of 771 from zero value for the classical case is due to the admixture of other 
fragments. 

Let us now consider more general characteristics of correlations by summing 77^ over a 
certain interval of i. In particular, after summing over all m fragments in the numerator 
and denominator of (|105| ) we have 

E E ({^pi) 2 ) A - (m) A 
V (m) = ^ . (109) 

Analogously, the summation over IMF's yields 

E E {(^vm)-)A - (M)a 
" (M) = — (M U (U0) 
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where = 6 and [12 = 40 were chosen. 

As follows from the results presented in Fig. 5, the above mentioned properties of the 
occupation number distributions survive in more global characteristics i](M) and r](m). 
Since the contribution of nucleons is dominant, therefore 77 (m) rs rjx for temperatures higher 
than about 2 — 3 MeV. Note that 77(77?,) goes down when T — > 0. Under this condition, the 
system of A identical particles exists as a unique fragment and r\A goes to its limiting value 
— 1 for all three statistics as discussed above. In other words, in the limit T — > it behaves 
like a closed-packed system ||28|| . It is natural that all correlation effects caused by quantum 



statistics are seen more clearly in the denser systems corresponding to smaller values of V. 

Rather similar correlation characteristics can be defined for the multiplicity distribution 
of a fragment i 

lt = A(N?) A -(N t ) A (111) 
where the dispersion A(N?)a is calculated according to Eq. (|89"D. Note that as compared to 



the definition of r\i 1 the factor (iVj) is dropped and the quantity 7^, which is the "shifted" 
multiplicity dispersion, may be considered as a shape distribution characteristics. 

The calculated 7« are displayed in Fig. 6. For the case of Poisson distribution, the 
dispersion (^) equals the mean multiplicity (Q) and therefore 7^ = 0. Such a situation 
is realised for a dominating nucleon component in the Boltzmann gas mixture where 71 ~ 
0. Visible differences appear for quantum cases. However, because in the whole mixture 
of different (non-identical) fragments the fraction of specific identical fragments is small, 
therefore all three statistics give similar results. As can be seen from Fig. 6 for i 7^ 1 and 
temperature just above the threshold, ji{T) is positive and has a spike-like structure showing 
that the fragment mass distribution is getting much wider at the temperature around the 
phase transition. At higher temperature, 7i(T) is negative, i.e., the distribution is sub- 



Poissonian |22j and approaches asymptotically the Poisson distribution. 



By suming over i one can define a more global characteristics for IMF's 



/'2 



— ^2 A(N?)a — (M)a (112) 



and for all produced fragments 



7 (m) = £ m?)A-{m) A . (113) 



i=l 



As is seen in Fig. 7, the above mentioned features are found in the 7(M) and 7(771) re- 
duced dispersions. The IMF distributions near the transition temperature have an enlarged 
multiplicity width and this width increases with increasing the freeze-out density of the 
fragmenting system. 



B. Thermodynamic quantities 

Let us now express the global thermodynamic averages characterizing the equation of 
state of nuclear matter through the partition functions Z A . The average energy of the system 
is 
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E = T 



2 d\nZ A 



dT 



v 



(114) 



After the substitution of the average energy (|38|) into expression for the ensemble averaging 
(|75|), we have 



A [A/i] 



i)A - E E £ a Z 



Z A i=1 l=1 



A-il 



(115) 



Here Eq. (|S3|) was used in getting the second equality and the function eTi was calculated 
according to Eq. (|55|) . Making use of Eq. (|84"D , the average energy in the system ( [L15| ) can 
be related to the multiplicity of fragments and to their average energies 



E = E el (N^ 



8=1 

x [A/i] 

Si = = X! ^ Za-u 

Ji 1=1 
_ [A/i] 

/i = E fa Z A-ij ■ 



(116) 
(117) 
(118) 



The evaluation of functions en and fu is considered in the Appendix B. 
The pressure in the multifragmenting system 



P=T 



d\nZ A 



dV 



dE 



{ vp^N^.^NA 
dV 



L/C to. 



pi 



dV 



(ypiiA 



(119) 



can be related to quantities averaged over single-particle momentum, if Eq. ([83]) is used for 
the mean occupation numbers 



i A [A/i] A 

p = y- E E p u z a-u = E Pi 

A i= i ; =1 7 = i 

-fj = === E ^A-i2 
/» (=1 

with the function P u calculated using Eq. fl5"9"|). 

The heat capacity at the constant volume can be found in the similar way 



(120) 
(121) 



dE 



^{N PI }N 1 ,...,N A ^ ~ ^{%}A^,...,A^>a 



J^2 



= P 2 E ^{vpiVp j) a £pi £pj ■ [ 
Averaging (|122 ) over momentum and using Eqs. ([S9|), (jB3|) and flS?D, we finally obtain 



122) 
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C v = (3 2 



-E 



1 



EE 



A—ik—jl 



+ 



E' 

i,l 



E il 



■>A-il 



(123) 



where e~a and e\ are defined by Eqs. fl55]) and (|63|), respectively. 

The classical limit of the derived thermodynamics relations can be found according to the 
procedure described above. The average energy in the QSM for the Boltzmann fragments is 
obtained from Eqs. (|116 ), ( 117 ) after the Boltzmann substitution fl24]) 



where e, is defined by 



(124) 



The pressure is given by Eqs. (|12(J|), (|121|) : 



E Pi w 



(125) 



i=i 



and Pi is defined by (ffl|). After applying the same procedure and using Eqs. (|8~9"D, (|101| ), 
( TO , (HH), the heat capacity (|23|) reduces to 



C* = (3 2 £ A B (NiNj) A el e] + (3 2 £ (e? - £" 2 ) (iV,) 



(126) 



i=l 



where ef is given by Eq. (|73|) . 

An example of thermodynamic characteristics calculated within the QSM of multifrag- 
mentation is presented in Fig. 8. There is a narrow region of temperatures where the average 
energy per nucleon in the system, E/A, is growing very quickly. This region clearly corre- 
sponds to a sharp maximum of the heat capacity Cy at a constant volume indicating the 
first order phase transition. The finite width of Cy reflects the finite size of the system under 
consideration. The position of the maximum of Cy depends noticeably on the fragmenting 
volume and is correlated with the position of maximum in IMF multiplicity (see Fig. 2). 
One should note that the quick increase in E/A as a function of T is unambiguously related 
to the "plateau" region of the inverse function T(E/A) which is called sometimes the caloric 
curve [] ||42|| . As follows from Fig. 8, the plateau region is more distinct for heavy systems. 

The volume dependence of the pressure is shown in Fig. 9. The presented family of 
isotherms clearly indicates the first order phase transition in the fragmenting system. Such 
a behavior in the QSM is caused solely by the Coulomb interaction of fragments. If the 
Coulomb potential ( [75] ) is neglected, then the isotherms follow the case of classical ideal 
gas (see r.h.s. in Fig. 9). One should stress also that the considered characteristics of the 
equation of state are insensitive to quantum statistics effects. 



1 To be exact, the caloric curve is the T(E*/A) dependence but the excitation energy per nucleon 
E* / A differs from E/A by a constant value which is equal to the binding energy of the fragmenting 
nucleus. 
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C. IMF multiplicity distribution 



In the canonical ensemble, a nuclear system of A nucleons may have A different species of 
fragments : 1 < i < A. A number of intermediate mass fragments M involving : \i\ < % < ^t 2 > 
is defined as 

M = Y,N t (127) 

where i is taken from the IMF range [i = (fi±, . . . , 1x2). So, a number of IMF's for the system 
of A nucleons may take the following values 

M = 0,1,...,[A/H . (128) 

The probability to find exactly M fragments in the mass range \x is 

Pa(M) = (129) 

where the partition function Z4 is given by Eq. ([H]) and the new partition function Z A (M), 
characterizing the probability of a system to be in a state with M, is defined as follows 

Z A {M)= J2 <*(£ Ni-M] n Q(N{) . (130) 

Let us derive the recurrence relations that are needed for calculating this new partition 
function Z A {M) . 

By the direct account of the restriction ( |127| ), let us generalize the mean occupation 
numbers of the A nucleon system to the case of M IMF's 

Ma(m) = —L-r J2 5 ( J2 n - m ) M Ni ,...,n a n QW (131) 

i 

where Q(Ni) is related to the variable xp t by Eq. @). By differentiating the logarithm of 
the partition function ( |130| ) with respect to the definition (§)), we get for the mean 

occupation numbers (|131|) : 

. \ / \ d\nZ A (M) , , 

(^) A (M) = Xpi j-^-^ . (132) 



The product of the partition functions Q(N{) for subsystems of identical fragments 
enters the definition of the partition function (|130|) . Each term of this product satisfies the 
recurrence relation fllP]), so we have 



x 



dQ(Ni 



^ (Tl) m 4, Q(Ni-l) , i = 1, 2, . . . , A . (133) 
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The direct differentiation of the partition function ( |130| ) with respect to xpi gives for all 
allowed values of M 



x. 



dZ A {M) 



pi 



dx 



pi 



E ME N k -M\ Up 



d\nQ(Ni 



dx 



n qw ■ ( i34 ) 



pi 



Substituting ( |133| ) into ( |134| ) and using the method employed above (see Eqs. 
to calculate (z^), one finds for the case M = 



x _ az A (M - o) = ^ ^ ZA ii(M = o) 

aa -p* 1=1 



X 



dZ A (M = 0) 



i £ fx 
i E \l 



By differentiating ( |134j ), we obtain in the case of M ^ and i E 



P)7 ( M\ min([A/i],M) 



<9x 



)-© ) 

(135) 
(136) 



(137) 



The mean number of fragments % for the A nucleon system with the fixed number M of 
IMF's is 



(N i ) A (M) = J2 Ma(M) 



(138) 



where the mean occupation numbers are determined by ( |131| ), ( |132| ). The mean fragment 
number (|138|) satisfies the restrictions (|36|) and (|127|) . By the substitution of Eqs. (|135|) , 
( |136| ) into ( |132[ ) and then into Eq. ( |138| ) we have for the mean fragment multiplicity in the 

case of M = 



(N t ) A (M = 0) 



[A/i] 



Z A (M = 0) 



E fuZ A _ u (M = 0) 



(Ni) A (M = 0) = 



i £ V 
i E ji . 



(139) 
(140) 



Here fu is calculated according to fl45|) . 

To find the average fragment multiplicity in the case of M ^ and i E fi, one inserts 
Eq. (|137| ) into (|132|) what after using (|138|) results in the following relation : 



(Ni) A (M) 



1 



mm([A/i},M) 

MM) E fuZ A -«(M-l) 



141) 



In accordance with ( |130| ), we have the following relations. For systems with the nucleon 
number j < \l\ and the IMF number M = 



Z j (M = 0) = Z j , i = 0, 1, . . . , ^ - 1 

For system with j < /ii • M and M^0 one has 



(142) 
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Z j (M) = , j = Q,l,...,n 1 -M-l . (143) 

Finally, for the "conditional" partition function Zj(M) in the case of M = 0, after the 
substitution of the mean fragment numbers ( |139| ), ( |I41| ) in the the conservation law of the 
total number of nucleons (0), we have the following recurrence equations 

i \JM 

Z j {M = 0) = - £ £ i /« ^(M = 0) > j = fJii,...,A . (144) 

3 i<£fi 1=1 

The recurrence equations for the case of any M ^ 0, (M = 1,2, ... , [j//ii]) are obtained 
similarly by making use of the restriction ( p.27| ) and Eq. ( |141| ) : 

, min{[j/i],M) 

Zj{M) = — E E fu Zj-u{M -I), j = f i 1 -M,...,A . (145) 

tG/i 1=1 

If one now divides Eqs. (|142|) and (|145|) by the partition function Zj and applies the prob- 
ability definition (|129|) , we get the recurrence relations for probability to find exactly M 
fragments among IMF's. If M — 



P,(M = 0) = 1 , j = 0,l,..., f i 1 -l (146) 

1 V/i] <v 

P j {M = 0) = - EE P j -ii(M = 0)if il ^± , j= f i l ,...,A . (147) 
■7 i^/i l=i L i 

If M ^ then 

P,(M) = , j = 0,l,...,fa-M-l (148) 

j min([j/i],M) „ 

p jM = mY, E Pj-uiM - 1) f u ^ , j = fa- M,...,A . (149) 

£GjU (=1 3 

In the Boltzmann limit, the IMF multiplicity distributions and related equations can be 
easily derived from the above derived quantum results by means of the substitution (p4]) . 
In particular, the equalities ( |142| ), ( |143| ) for the partition function are not changed and 
Eqs. (|144j) and ( |145|) reduce to the following equations for M = and M^O, respectively : 

Z?{M = 0)= - E ku k Zf_ k (M = 0) , j = fa,...,A (150) 

Zf{M)= ±- " k Z?_ k (M - 1) , j = fa- M,...,A . (151) 

Similarly, the recurrence equations for the probability to find M IMF's in the classical model 
of multifragmentation result from Eqs. ( |146p — (|I4"9"D . If M = 

Pj(M = 0) = l , j = 0,1,... ,1^-1 (152) 

p.( M = 0)= - E Pj-k(M = 0) k u k , j = Hi,...,A . (153) 



J 
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If M^O 



Pj{M)=Q , j = 0,l,...,fi X -M-l (154) 

PjW) = i £ P/-*(M - 1) ^ , j = f i 1 >M,...,A . (155) 

One should note that the multiplicity distributions of IMF's within the Boltzmann statistics 
have been recently studied in Refs. j|3,H|. Their equation coincides with our Eq. (|155 ) 



but Eqs. (|152|) , ( |153| ) for M = have not been found and the normalization condition 
£ M P(M) = 1 were used to get the full distribution in |43],|44j]. 

The multiplicity distributions of IMF's (/ii = 6, /i 2 = 40) calculated within the QSM of 
multifragmentation are presented in Fig. 10. It is of great interest that the M-distribution 
sharply changes its shape in the vicinity of the phase transition temperature defined by the 
maximum of the heat capacity. These results are in full agreement with the discussion of the 
reduced dispersions in Figs. 6, 7 and reveal the different IMF distributions in different phases 
as it would be expected for the phase transition. Details of the multiplicity distributions 
depend on the mass number A and the size of the system, being weakly sensitive to quantum 
statistics of fragments. Some local irregularities and the appearance of two maxima in the 
transition curve of IMF multiplicity distributions are caused by the use of experimental 
values for the fragment binding energy. If the Coulomb fragment interaction is neglected 
and a smoothed approximation is used for the binding energies, the QSM with the Boltzmann 
statistics reproduces exactly the numerical results in [I3,31|. 



D. Total multiplicity distribution 

The total fragment multiplicity m associated with the multifragmentation of a system 
of the size A is defined similarly to ( |127| ) 

m = X>< (156) 

except that the region fi is now : 

H = (1,2,...,A) (157) 

i.e., the events with m = are absent due to the nucleon number conservation. The proba- 
bility to find m fragments is defined by the same Eq. (125) with the conditional probability 
( |130|) . From this follows immediately that in the case of m = the only non-zero partition 
function is that for j = : 

Z j (m) = 5 jfi , j = 0,1,..., A . (158) 

For the partition functions with j < m, the following relations are valid 

Z j {m)=0 , m = 1,2, ...,A, j = 0, 1, . . . , m - 1 . (159) 

In this case, the recurrence equations ( |133| ) and the derivative of the partition function ( |134| ) 
still hold but the region \i is defined by Eq. ( |157 ). Similar transformation of Eq. ( |134j ) with 
using (|133|) gives for the particular value m = 1 
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X 



pi 



X 



d Z A (m) 

dxpi 
dZ A (m) 



pi 



dx 



xpi Z A „i(m - 1) 



i = l,2, 



i = A . 



A - I 



(160) 
(161) 



pi 



For derivatives ( |134j ) with m = 2, 3, . . . , A we have 



.r 



OXpi 

dZ A {m) 



m'm([A/i],m) 

]T ( T l) i+1 4, ^_«(m - , z = l,2,...,A-m + l 



pi 



dx 







i = A- m + 2,...,A . 



(162) 
(163) 



pi 



According to the definition, the mean fragment number of the A nucleon system with the 
total multiplicity m, for which two restrictions (^) and ( |156| ) are fulfilled, is 



(iVi)A(m) = £ (%)a( 



m 



(164) 



By substituting ( |160| ) and ( |161| ) into Eq. ( |132j ), we get for the mean fragment number ( |164| ) 
in the case of m = 1 : 



(N l ) A (m) = , 
(JVi) A (m) = /a 



[m — 1) 



Za{™) 



i = 1,2,..., A- 1 



(165) 
(166) 



with /jx calculated according to fl45|) . This quantity for the case of m = 2, 3,..., A is 
obtained by combining ( |162| ), ( |163| ) and ( |l32p 



(iVi>A(m) 



1 



min( [A/i],m) 

E fuZ A _ a {m-l) , z = l,2,. ..,A-m + l (167) 
(JV i ) i4 (m) = , 2 = A-m + 2,...,A (168) 



with defined again by (|45|) . 

As usually, the recurrence relations are derived from the restriction ([36|) to the j nucleon 
system. From Eqs. ( pL65| ) , ( |166| ) we have for m = 1 



Zj(m) 



fji , j = l,...,A (169) 

and similarly for m = 2, 3, . . . , A, by using Eqs. ( |167|) and (|168|) , we have 



Zj{m) 



j—m+l min([j/i],m) 

7 X) ^2 i fu Zj-uim - /) 



j = m, . . . , A . 



(170) 



Eqs. ( |158j ), (|159|) and ( |169| ), (|170|) form a complete set which is sufficient for calculating the 
partition function Z A (m) and therefore the multiplicity distributions P(m). 

The limiting case of the Boltzmann statistics for the total multiplicity distributions can 
be obtained from the appropriate quantum equations by means of the substitution 
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For the system of classical fragments, Eqs. (|158|) , ( |159|) , ( |169|) are not changed and Eq. fl2"4|) 
results in the following recurrence relations in the case of m = 2, 3, . . . , A : 

^ j-m+l 

Zf{m)= - J2 iUiZf^m-l) , j = m,...,A. (171) 

J i=i 

The total multiplicity distributions of fragments within the QSM are exemplified in 
Fig. 11. General behavior of m-distributions is close to those noted before while discussing 
the IMF distributions (cp. to Fig. 10), however all curves became more smooth. The reason 
of this behavior is a dominance of individual nucleons what results also in somewhat stronger 
influence of the quantum statistics on the multiplicity distributions. 

V. CONCLUSIONS 

The model proposed in this work realizes for the first time the quantum statistical ap- 
proach to the multifragmentation of excited finite nuclei in the framework of the canonical 
ensemble method. The mathematical basis of the model is the recurrence equation tech- 
nique which allows, within the usually accepted physics assumptions, to solve a cumbersome 
problem of fragment partitioning exactly without involving complicated and time-consuming 
Monte Carlo methods. The QSM opens a possibility for calculating in the same framework 
various microscopic characteristics of occupation numbers, global thermodynamic variables 
specifying the equation of state, and different observables allowing for a comparison with 
the experiment. The exactly solvable technique makes such calculations transparent to the 
physics assumptions used. The QSM includes the limiting case of the classical Boltzmann 
statistics and in this case the QSM reduces to a class of solvable models developed ear- 

her fgrrjjBjg. 

In principle, effects of quantum statistics in the QSM are seen on the microscopic level in 
quantities like (ui), (z/ 4 2 ), especially at temperatures corresponding to the maximal fragment 
yield. The strength of these effects is proportional to a number of available identical clusters 
of a given size. Thus, the effect is strongest for nucleons and is noticeably weaker for the 
IMF's. It is of interest that large uncertainties in the choice of the freeze-out volume do not 
kill these quantum statistical effects in the correlation characteristics fy. 

Global thermodynamic variables, which are insensitive to quantum statistics, exhibit 
nevertheless a rather sharp phase transition of the first order in the temperature region w 6 — 
8 MeV depending on the size of the system and on the value of the freeze-out volume V/Vq. 
In the QSM, the observed Van der Waals behavior is caused by the competition between 
attractive forces, included in the fragment binding energies, and the Coulomb interaction 
energy. It is noteworthy that at not so high temperatures, the value of V/Vq = 3, being 
the standard value for the SMM calculations of the primary fragment partition |J, turns 
out to be in the thermodynamically unstable region : dP/dV < 0. This implies that before 
getting an equilibrium solution in the model analysis, the Maxwell construction should be 
applied to the isotherms in this unstable region. The use of these Maxwellian isotherms is 
important for extracting true values of the freeze-out temperature of a system as well as for 
detailed studies of the caloric curve. 

Fragment multiplicity is a very interesting observable, even though the quantum statis- 
tical effects are practically absent in the averaged quantities. The mean number (iVj) of 
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IMF's clearly demonstrates the "rise and fall" behavior as a function of temperature. Frag- 
ment multiplicity dispersion has a maximum near the transition temperature corresponding 
to the maximum of heat capacity Cy as a sign of the phase transition in finite systems to 
be dependent on the system size, freeze-out volume and fragment binding energies. These 
features are remarkably well seen in the change of the shape of the IMF multiplicity distri- 
bution in the transition region. The measurement of the HBT effect for IMF's which would 
be crucial for determining the size of a fragmenting system, seems to be a hard task due 
to a comparatively small number of identical particles. However, the model studies in this 
paper show where this effect is expected to be maximal. In addition, when searching for the 
HBT effect a more detailed quantity, the correlated particle spectrum (N^p) Ni(j?)) rather 
than the mean square multiplicity (Nf) should be investigated [p2| . In any case, this issue 
deserves a special discussion. Making use of protons instead of the IMF's unfortunately does 
not help much because in a real experiment this component may be strongly enriched by 
non-equilibrium particles. 

As a final remark, the model developed in this paper takes into account the mass number 
conservation while the conservation of the charge number may be quite important for the 
yield of specific isotopes or even for the equation of state. In particular, it was shown that 
the liquid-gas transition in the asymmetric nuclear matter is of the second order rather than 
of the first order as one would expect for neutron-proton symmetric matter [25]. This process 



is of particular interest because experimental data for Au + C at 1 GeV per nucleon [46 



shows a much smoother phase transition than the first observation in the Au + Au collision 



at 600 MeV j4T[. The recurrence equation technique allows one to consider simultaneously 



both mass and charge number conservation. Such a generalization has been done in Refs. 



[47,48] but only for the mean multiplicity of individual fragments within the Boltzmann 
statistics. Similar extension for calculating the multiplicity distributions of the identified 
IMF's would open a new way to study critical exponents, intermittence phenomenon and so 
on. This work is in progress now. 
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APPENDIX A: A USEFUL RECURRENCE RELATION 



As an auxiliary step, let us define the partition function Q n ^ for a system of A" Bose 
(Fermi) identical particles with the fixed number rip of particles in a given single particle 

state p 

Qn = E S ( E - N ) S far - "?) II x 7 ( Al ) 
{»?} \f J 9 

where the conservation laws are accounted for by the Kronecker symbol and the allowed 
occupation numbers up are given by (Ej). Then, the total partition function for A" 
identical particles can be expressed via auxiliary one (|A1[) as follows 

{l},{iV} 

Qtv= E (A2) 

rtp=0 

where the upper limit in the sum is equal to AT or 1 for bosons or fermions, respectively 
The sum in (|A1|) with respect to up is taken by means of delta functions 

q7= e * [ e »? + np-A x y n # 

= E 8 ( E ^ ~(N- np)\ * far) II # • ( A3 ) 
0y} Vp* / P' 

Comparing the r.h.s of ( |A5| ) to the definition (|A1[ ), one can see that here one can select out 
the partition function for (AT — np) particles in the state p = 

Q7 = #Ql£l ■ (A4) 



This will result in the recurrence relations for the partition function ( M] ) which are valid 



for both Bose and Fermi particles. Let us start with the boson case. After the substitution 
of ( |A4| ) into ( |A2| ), we have for the total partition function 



N 

Qn = Q7 =0 + E # • (A5) 



Changing the summation indices in the second term on the r.h.s. of (|A5|) and using the 
recurrence relations (|A4j) , we obtain 



N-l 

Qn = Q H n~ + xp E Qn-i = Qn 3 ' + x p Qn-i (A6) 

np~=0 

where the partition function for (N — 1) particles appears explicitly in the second term. 
After differentiating ( |A6|) with respect to xp and taking into account that this derivative of 
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the partition function Q"^ equals zero, we arrive at the following differential recurrence 
equations: 



p 



x p Qat_i + — . (A7) 



The mean occupation numbers can be expressed through the derivative of the logarithm of 
the partition function Q N with respect to xp (see Eq. fllOD ). So, Eq. ([AT]) for Bose particles 
reduces to the recurrence equations 

(rif) N Q N = xp Q N -i [1 + {np) N -i\ (A8) 

where Q = 1 and (np) = according to Eqs. (||) and (|7|), respectively. 

Let us proceed now to the system of fermions. Plugging ( [A"4] ) in ( |A2|) , we have for the 
partition functions of two systems with N and (N — 1) particles 

Qn = QT° + *pQTi (A9) 

Qn-i = QTi+^pQTi-i • (A10) 

To get the differential equations for the system of fermions, let us differentiate Eqs. (|A9|) 
and ( |A10| ) with respect to xp, multiply the first equation by xp and second one by x| and 
sum these two : 

dQ N _ , <9Qjy-i . , 

Xp> ^ - XpUi N -l Xp . (Allj 

Here we took into consideration that derivatives of Q^T ° and Q^-i with respect to xp 
equal zero. Combining ([H^) and ( |A11|) , we obtain the final recurrence equations for the 
Fermi system : 

(np) N Q N = xpQ N _! [1 - {np) N -i\ . (A12) 

The same initial condition as in the boson case are used here. One should note that Eqs. 
( A"8|) and ( |A12|) for Bose and Fermi statistics differ only by the sign in front of the second 



term. 



APPENDIX B: MOMENTUM RESUMATION 



Using (|T1D, let us present the sums in (| 

p 



and 



fa = (Tiy +i e 



Pa 



l+l c -PlEi 



in more convenient form 



(Bl) 



, y +1 m y> ]?_ 



pip 2 

2m, 



+ Ei fa 



(zpl)W e"^ ]T 
p 

3 V ~ ;<l V 



^2 \2 

' 1 e 2™, + 2£i^ 



2m,- 



2 

+ 317, 



(B2) 

(B3) 
(B4) 
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pip 2 


= ^ 6 " 




p 


ri 


f pip 2 


_ (Ap) 2 


2rrii 


2mj 



The summation over momentum in the above Eqs. (B1|)-([B4|) can be written as 

2m; (B5) 
f3l(Ap) 2 ->2 

y — e~2m i = ^LL V n 2 e~ ^ U (B6) 
2m; - 

£ (4) e ~ 2m -y-£-) ?* e ~ 2m ' • (B7) 

where n = (n^,n^,u^) and each nS a ' runs over all integer numbers (see Eq. (P). It is 
worthwhile to proceed to the momentum integration in (|B5|) -( [B7|) , though this is not valid 
in the whole region of temperature T and volume V. The passage to the integration over 
momentum 



is possible under the condition 



In the opposite case, when 



(Apf 



1. (B9) 
2m,-T ~ v ; 



> 1 (BIO) 



this passage can be erroneous and the sums (|B5|)— flB7| ) should be calculated numerically. If 



the condition (fB9| ) is satisfied, we have for integrals of (p5|)-(p7|) : 



l y 

^3 (bid 



v 



y P_ e 2m, = - T -A- (B12) 
^ 2m,- 2 / 5 /2 A? 

p 



D 2 3 „ 1 V 



S(4 K^T^I < B13 > 



where the thermodynamic wavelength for a fragment made from i nucleons is 

, /2^ 2 \ 1/2 , . 

Ai = — — , rrii = im N (B14) 

\rriiT J 

and 77i 7v stands here for the nucleon mass. After the substitution of ( Bll|) -( |Bi3|) into ( Bl|) - 
£3[), we obtain 
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1 V 



h = (=F1)' +1 «-<"* W jf 



-TT 1 + Ei) fa 



Z_ T 2 r 2 + 3 £. T r l + ^2 ) fa 



(B15) 
(B16) 
(B17) 



For calculating (|B5|) (|E7|) in the case ( |B10| ), let us present these sums in the form more 
suitable for computing. Let us introduce the function depending on the integer vector n 



zu(<r)= J2 n ° e ~ 2mi n > a = 0,2,4 . 

n=— oo 

In this notation, the sums (p5|)-(p7|) are reduced to the following expressions : 

4(0) 

(Ap) 2 



(B18) 



E 



P 



-HO 

2m,- 



e = 3^(2)^(0) 



2m, ; 



E 



' \ 2 r 

e 2m s = 3^(0) [^(4) ^(0) + 2 4(2) 



(Apr 

4m? 



(B19) 
(B20) 
(B21) 



To calculate them, it is sufficient to know functions zu(0),zu(2) and ^(4) in the limited 
range of values n, because exponents go fast to zero with increasing the absolute value of n. 
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Figure captures 



Fig. 1. The mean occupation numbers (ypi) of fragments of mass i for A = 200 system 
at temperatures of 10 and 15 MeV. The dashed, dotted and continuous curves correspond 
to the QSM result obtained under the assumption that all fragments are bosons, fermions 
or Boltzmann particles, respectively. The system volume is V/V = 3 and Ap = 90 MeV/c. 

Fig. 2. The temperature dependence of the mean multiplicity (Nj) for fragment of mass 
% in the fragmenting system of size A = 200. For further details, see the caption of Fig. 1. 
The results are obtained for the two values of the system volume. 

Fig. 3. The temperature dependence of the total mean multiplicity (m) and IMF multi- 
plicity (M) for systems with A = 100 and 200. For further details see the caption of Fig. 2. 

Fig. 4. The temperature dependence of the occupation number correlation rji of the 
fragment species i for the system with A = 200 at two different values of parameter V/Vq. 
For further details see the caption of Fig. 2. 

Fig. 5. The temperature dependence of the occupation-number correlation parameter 
summed over all fragments r](m) and IMF's r/(M) for A = 200. For further details see the 
caption of Fig. 4. 

Fig. 6. The temperature dependence of the shifted multiplicity dispersion 7, of the frag- 
ment species i for the system with A = 200 at two values of the freeze-out volume parameter 
V/Vq. For further details see the caption of Fig. 2. 

Fig. 7. The temperature dependence of the shifted multiplicity dispersion summed over 
all fragments 7(m) and IMF's 7(M) for the system with A = 200. For further details see 
the caption of Fig. 4. 

Fig. 8. The heat capacity and the total energy vs. temperature for systems with A = 200 
(solid line) and 100 (dashed) nucleons at two different values of the volume parameter V/Vq. 
Dotted curves are calculated for A = 200 system neglecting the Coulomb interaction energy. 
The results for Bose, Fermi and Boltzmann statistics coincide with each other with the ac- 
curacy comparable to the line thickness. 

Fig. 9. The pressure as a function of the volume at the constant temperature T for the 
systems size A = 100 and 200 with (l.h.s.) and without (r.h.s.) taking into account the 
Coulomb interaction energy. The results for Bose, Fermi and Boltzmann statistics coincide 
with each other with the accuracy comparable to the line thickness. 

Fig. 10. The IMF multiplicity distributions for systems with A = 200 and 100 nucleons 
at the given temperature T for two different values of the volume parameter : V/Vq = 6 
(lower panels) and V/Vq = 3 (upper panels). Open squares, open circles and filled triangles 
are the QSM result obtained under the assumption that all fragments are bosons, fermions 
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or Boltzmann particles, respectively. 

Fig. 11. The total fragment multiplicity distributions. For further details see the caption 
of Fig. 10. 
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